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ABSTRACT 

The properties of the chiral phase transition at finite temperature and 
chemical potential are investigated within a nonlocal covariant extension of 
the Nambu-Jona-Lasinio model based on a separable quark-quark interaction. 
We consider both the situation in which the Minkowski quark propagator has 
poles at real energies and the case where only complex poles appear. In the 
literature, the latter has been proposed as a realization of confinement. In 
both cases, the behaviour of the physical quantities as functions of T and 
jj, is found to be quite similar. In particular, for low values of T the chiral 
transition is always of first order and, for finite quark masses, at certain "end 
point" the transition turns into a smooth crossover. In the chiral limit, this 
"end point" becomes a "tri critical" point. Our predictions for the position of 
these points are similar, although somewhat smaller, than previous estimates. 
Finally, the relation between the deconfining transition and chiral restoration 
is also discussed. 
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The behaviour of hot dense hadronic matter and its transition to a plasma of quarks 
and gluons has received considerable attention in recent years. To a great extent this is 
motivated by the advent of facilities like e.g. RHIC at Brookhaven which are expected 
to provide some empirical information about such transition. The interest in this topic 
has been further increased by the recent suggestions that the QCD phase diagram could 
be richer than previously expected (see Ref.[|IJ for some recent review articles). Due to 
the well known difficulties to deal directly with QCD, different models have been used 
to study this sort of problems. Among them the Nambu-Jona-Lasinio model is one 
of the most popular. In this model the quark fields interact via local four point vertices 
which are subject to chiral symmetry. If such interaction is strong enough the chiral 
symmetry is spontaneously broken and pseudoscalar Goldstone bosons appear. It has 
been shown by many authors that when the temperature and/or density increase, the 
chiral symmetry is restored [Q]. For zero chemical potential and finite temperature this 
transition is found to be a smooth one. However, whether for finite chemical potential 
and zero temperature the transition is of first order or not is highly dependent on the 
parameters of the model and on the approximations made. For example, in the SU(2) 
version of the NJL model a Hartree-Fock treatment leads to a first order transition, while 
the Hartree treatment as well as the SU(3) model calculations seem to favor a second 
order one||. On the other hand, several recent calculations performed within various 
other models^, ||, [J clearly suggest that in QCD this transition should be of first order. 
Interestingly, the model used in Ref.|| can be understood as a nonlocal generalization 
of the NJL. However, the corresponding four point vertex is local in time and, therefore, 
not covariant. Some covariant nonlocal extensions of the NJL model have been studied 
in the last few years @. Nonlocality arises naturally in the context of several of the most 
successful approaches to low-energy quark dynamics as, for example, the instanton liquid 
model and the Schwinger-Dyson resummation techniques [[|]. It has been also argued 
that nonlocal covariant extensions of the NJL model have several advantages over the 
local scheme. Namely, nonlocal interactions regularize the model in such a way that 
anomalies are preserved [fT0f| and charges properly quantized, the effective interaction is 
finite to all orders in the loop expansion and therefore there is not need to introduce extra 
cut-offs, soft regulators such as Gaussian functions lead to small NLO corrections [[□]], etc. 
In addition, it has been shown [|I2}| that a proper choice of the nonlocal regulator and the 
model parameters can lead to some form of quark confinement, in the sense of a quark 
propagator without poles at real energies. Recently, the behaviour of this kind of models 



at finite temperature has been investigated [13]. In this work we extend such studies to 
finite temperature and chemical potential. 

We consider a nonlocal extension of the SU(2) NJL model defined by the effective 
action 

S = d A x tp(x) (ift — m c ) ip(x) + / d A xi...d 4 X4 : V(x±, X2, X3, X4) 
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x (i)(xi)i)(x3)i)(x 2 )i)(x 4 ) + i){xi)i^ h T a ^{x^{x 2 )i^T a ^{x A )j , (1) 

where m c is the (small) current quark mass responsible for the explicit chiral symmetry 
breaking. The interaction kernel in Euclidean momentum space is given by 

V{q x , q 2 , q 3 , 94) = G r(qf)r(ql)r{ql)r(ql) 5{q x + q 2 -q 3 - qi) , (2) 

where r(q 2 ) is a regulator normalized in such a way that r(0) = 1. Some general forms 
for this regulator like Lorentzian or Gaussian functions have been used in the literature. 
A particular form is given e.g. in the case of instanton liquid models. 

Like in the local version of the NJL model, the chiral symmetry is spontaneously 
broken in this nonlocal scheme for large enough values of the coupling G. In the Hartree 
approximation the self-energy S(g 2 ) at vanishing temperature and chemical potential is 
given by 

£(g 2 ) =m c + (£(0)-m c )r 2 (g 2 ), (3) 
where the zero-momentum self-energy S(0) is a solution of the gap equation 

^ wo) - m c) = ! A ^iMiBWra . (4) 
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In general, the quark propagator might have a rather complicate structure of poles and 
cuts in the complex plane. In what follows we will assume that the regulator is such that 
it only has an arbitrary but numerable set of poles. As already mentioned, the absence 
of purely imaginary poles in the Euclidean quark propagator might be interpreted as a 
realization of confinement HT2| . In that case quartets of poles located at a p = R p ± i I p , 
a p = —R p ± i I p appear. On the other hand, if purely imaginary poles exist they will 
show up as doublets a p = ±z I p . It is clear that the number and position of the poles 
depend on the details of the regulator. For example, if we assume it to be a step function 
as in the standard NJL model only two purely imaginary poles at ±i M appear, with M 
being the dynamical quark mass. For a Gaussian interaction, three different situations 
might occur. For values of S(0) below a certain critical value £(0) cr ;t two pairs of purely 
imaginary simple poles and an infinite set of quartets of complex simple poles appear. 
At S(0) = S(0) cr j t , the two pairs of purely imaginary simple poles turn into a doublet 
of double poles with R p = 0, while for E(0) > E(0) cr j t only an infinite set of quartets of 
complex simple poles is obtained. For the Lorentzian interactions there is also a critical 
value above which purely imaginary poles cease to exist. However, for this family of 
regulators the total number of poles is always finite. 

To introduce finite temperature and chemical potential we follow the imaginary time 
formalism. Thus, we replace the fourth component of the Euclidean quark momentum 
by oj n — ifi, where u n = {2n + 1)ttT are the discrete Matsubara frequencies and \i is the 
chemical potential. In what follows we will assume that the temperature and chemical 
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potential dependences enter only through those quantities that carry a (^-dependence at 
T = fi = 0. That is to say, we will consider that the model parameters G and m c , as well 
as the shape of the regulator, do not change with T or \x. Performing this replacement in 
the gap equation Eq.(|J) we get 



2 7T 4 

GN n 



„ n=oo 

:S(0)-m c ) = 27rT J d 3 q £ F(q iu>n ) 



where 



and q 2 
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ifi) 2 + q 2 . The sum over n can be expressed in terms of a sum over 
the poles of J-(q z ), that we denote z p , by introducing the auxiliary function f(z) = 
1/(1 + exp(z/T)) and using the standard techniques described e.g. in Ref . |fL4|| . It is easy 
to see that the poles of J-(q z ) are closely related to those of the quark propagator a p . For 
the pair a p = ±(R P + il p ) the associated values are z p = ±e p — /x =F iRpIp/t P , while for 



±(R P — il p ) one has z p = ±e p — \i ± iR p I p /e p , where e p is defined as 
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Assuming that all the poles are simplef] and performing explicitly the sum over all the 
poles in a given multiplet we get 
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where ResftzF^z); a p ) stands for the residue of the function 1zT{z) evaluated at 
and the prime in the sum indicates that it runs over all the poles a p = R p + il p with 
R p > and I p > 0. The coefficient 7 P is defined as 7 P = 1/2 for R p = and j p = 1 
otherwise. The generalized occupation numbers n± are given by 

-i 
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As customary, in writing Eq.(|8|) we have isolated a term which has the same form as the 
T = fi = expression. In this way, all the T and \i dependent contributions remain finite. 
Replacing Eq. (||) in the right hand side of Eq. (|5|) and using 

nu)r 2 (u) 



Res[2,2 J-(z); a p ] 



d u mu) 



(10) 



lr This assumption is made just to keep our expressions into a simple form. The generalization to poles 
of arbitrary multiplicity is rather straightforward. 
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we obtain the final form of the gap equation at finite temperature and chemical potential. 
As we see the dependence on these quantities is completely fixed by the pole structure of 
the quark propagator. 

The other quantities which are of interest to understand the characteristics of the 
chiral phase transition are the quark condensate (qq) and the quark density (q'q) for each 
flavour. At T = fi = 0, the condensate is given by 

d 4 q 



(qq) 



Tr 



+ £(<? 2 
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Following similar steps as before, the corresponding result for finite T and \x can be cast 
into the form 



{qq) 



4vr 4 



£(g 2 
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(12) 



Away from the chiral limit this expression turns out to be divergent. Following the 
standard procedure, we regularize it by subtracting the value obtained in the absence of 
interactions. 

To determine the quark density one has to be more careful due to the presence 
of nonlocal interactions. In fact, they imply the existence of extra contributions to the 
Noether currents. The proper expression for the quark density in the Hartree approxima- 
tion reads 

d 4 q 



(qiq) = -i 



(27T)' 



Tr 



(13) 



It is seen here that, for any quark flavour, the residue of the pole given by the dressed quark 
propagator is equal to one. As shown in Ref. |]15||, this leads to the correct normalization 
of the baryon number, independently of the shape of the regulator. As above, the result 



in Eq.(13) can be now extended to finite temperature and chemical potential by replacing 
g 4 — > uj n — i \i and summing over all Matsubara frequencies cu n . In this case we obtain 
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Having introduced the formalism needed to extend the model to finite temperature 
and chemical potential we turn now to our numerical calculations. In this work we take 
the nonlocal regulator to be of the Gaussian form 



r(q 2 



exp 



2A 2 



(15) 



and consider two sets of values for the parameters of the model. Set I corresponds to 
G = 50 GeV" 2 , m c = 10.5 MeV and A = 627 MeV, while for Set II the respective values 
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are G = 30 GeV -2 , m c = 7.7 MeV and A = 760 MeV. Both sets of parameters lead to the 
physical values of the pion mass and decay constant. For Set I the calculated value of the 
chiral quark condensate at zero temperature and chemical potential is —(200 MeV) 3 while 
for Set II it is —(220 MeV) 3 . These values are similar in size to those determined from 
lattice gauge theory or QCD sum rules. The corresponding results for the self-energy at 
zero momentum are S(0) = 350 MeV for Set I and S(0) = 300 MeV for Set II. Using the 
explicit expression of T,(0) crit for the Gaussian interaction, 



S(0) cri4 = m c + i (\fm 2 c + 2A 2 - m}j exp 



ml + 2A 2 + m, 



4 A 2 



(16) 



it is easy to check that Set I corresponds to a situation in which there are no purely 
imaginary poles of the Euclidean quark propagator and Set II to the case in which there 
are two pairs of them. Therefore, we will also refer to Set I as the confining set and to 
Set II as the non-confining one. 

The behaviour of the zero- momentum self-energy S(0), the chiral quark condensate 
(qq) and the quark density p = (q^q) as functions of the chemical potential for some values 
of the temperature is shown in Fig. 1. The left and right panels in the figure correspond 
to the results for Set I and Set II, respectively. In both cases we observe the existence of 
some kind of phase transition at (or around) a given value of the chemical potential which 
depends on the temperature. To obtain our results, we have included in the sums over 
a p appearing in Eqs.(|3|), flI2|) and flI3|) the first few poles of the quark propagator. We 
have checked, however, that for the range of values of T and \x covered in our calculations, 
the convergence is so fast that already the first pole gives almost 100% of the full result. 
Thus, the behaviour of relevant physical quantities up to (and somewhat above) the phase 
transition is basically dominated by the first pole of the quark propagator. 

We observe in Fig. 1 that at T = there is a first order phase transition for both 
the confining and the non-confining sets of parameters. As the temperature increases, the 
value of the chemical potential at which the transition shows up decreases. Finally, above 
a certain value of the temperature the first order phase transition does not longer exist 
and, instead, there is a smooth crossover. This phenomenon is clearly shown in the right 
panel of Fig. 2, where we display the critical temperature at which the phase transition 
occurs as a function of the chemical potential. The point at which the first order phase 
transition ceases to exist is usually called "end point" . In the chiral limit the latter turns 
into the so-called "tricritical point", which is the point at which the second order phase 
transition expected to happen in QCD with two massless quarks becomes a first order 
one. In fact, this is also what happens within the present model in the chiral limit, as it 
is shown in the left panel of Fig. 2. Some predictions 0, [f| about both the position of this 
point and its possible experimental signatures |16| exist in the literature. In our case the 



'tricritical point" is located at (T P ,/i P ) =(70 MeV, 130 MeV) for Set I and (70 MeV, 140 
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MeV) for Set II, while the "end points" are placed at (T e ,[me) =(70 MeV, 180 MeV) and 
(55 MeV, 210 MeV), respectively. As we see the predicted values are very similar for both 
sets of parameters and slightly smaller than the values in Refs.|5|, §], Tp ~ 100 MeV and 
/ip fa 200 — 230 MeV. In this sense, we should remark that our model predicts a critical 
temperature at /x = of about 100 MeV, somewhat below the values obtained in modern 
lattice simulations ||17|| which suggest T c « 140 — 190 MeV. In any case, our calculation 
seems to indicate that fip might be smaller than previously expected even in the absence 
of strangeness degrees of freedom. 

It is interesting to discuss in detail the situation concerning the confining set. In this 
case we can find, for each temperature, the chemical potential /x^ at which confinement 
is lost. Following the proposal of Ref.[|12[], this corresponds to the point for which the 
self-energy at zero momentum reaches £(0) cr jt (c.f. Eq.(|l6|)). Using the values of m c and 
A corresponding to Set I we get S(0) cr .j i = 267 MeV. For low temperatures, /x^ coincides 
with the the chemical potential at which the chiral phase transition takes place. However, 
for a temperature close enough to that of the "end point" , /x^ starts to be slightly smaller 
than the value of /x that corresponds to the chiral restoration. Above Tg it is difficult to 
make an accurate comparison since, for finite quark masses, the chiral restoration proceeds 
through a smooth crossover. However, we can still study the situation in the chiral limit. 
In this case we find that, in the region where the chiral transition is of second order, 
deconfinement always occurs, for fixed T, at a lower value of /x than the chiral transition. 
The corresponding critical line is indicated by a dashed line in the left panel of Fig. 2. In 
any case, as we can see in this figure, the departure of the line of chiral restoration from 
that of deconfinement is in general not too large. This indicates that within the present 
model both transitions tend to happen at, approximately, the same point. 

In conclusion, in this work we have investigated a nonlocal covariant extension of 
the NJL model at finite temperature and chemical potential. We have assumed that the 
corresponding four point vertex is separable in momentum space and that the regulator 
leads to an Euclidean quark propagator with an arbitrary (but discrete) number of poles 
in the complex plane. We have shown that the dependence on /i and T can be expressed 
in terms of these poles and the corresponding residues. We have studied in detail the 
particular case of a Gaussian regulator, considering two different situations. In the first 
case, the Euclidean quark propagator does not have any purely imaginary poles, which 
can be understood as a way of confinement. In the second case, this type of poles exists, 
thus quarks would be not confined. We have found that in both cases the behaviour of 
the physical quantities as a function of T and /x is quite similar. In particular, for low 
values of T the chiral transition is always of first order and, for finite quark masses, at 
certain "end point" the transition turns into a smooth crossover. In the chiral limit, this 
"end point" becomes a "tricritical" point. Our predictions for the position of these points 
are similar, although somewhat smaller, than previous estimates obtained in alternative 
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models^ |J. It is clear that several extensions of this work are possible. Firstly, it would 
be interesting to include interactions in the vector and axial-vector meson channels, in 
order to investigate meson properties close to the phase transition. Secondly, to study 
the possible existence of color superconductivity in this sort of models, interactions in the 
quark-quark channels should be incorporated. Finally, the extension to the SU(3) sector 
would allow to understand the role of the strangeness in the characteristics of the chiral 
phase transition as described by nonlocal covariant models. We hope to report on these 
issues in forthcoming publications. 
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Figure 1: Behaviour of the self- energy, the chiral condensate and the quark density as 
functions of the chemical potential for three representative values of the temperature. Full 
line corresponds to T=0, dashed line to T = 50 MeV and dotted line to T = 100 MeV. 
The left panels display the results for Set I and the right panels those for Set II. The quark 
density p is given with respect to nuclear matter density p = 1.3 x 10 6 MeV 3 . 
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Figure 2: Critical temperatures as a function of the chemical potential. The left panel 
corresponds to the chiral limit and the right panel to the case of finite quark masses. In 
the left panel, the dotted lines correspond to the region of the second order phase transition, 
while the full lines in both panels indicate the region where a first order transition occurs. 
The dashed line in the left panel indicates the critical line for deconfinement corresponding 
to Set I. For chemical potentials somewhat larger than /ip this line coincides with that of 
the chiral restoration. 



10 



